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In this paper we examine the use of an adiabatic quantum data transfer protocol to build a 
universal quantum computer. Single qubit gates are realized by using a bus protocol to transfer 
qubits of information down a spin chain with a unitary twist. This twist arises from altered couplings 
on the chain corresponding to unitary rotations performed on one region of the chain. We show 
how a controlled NOT gate can be realized by using a control qubit with Ising type coupling. The 
method discussed here can be extended to non-adiabatic quantum bus protocols. We also examine 
the potential of realizing such a quantum computer by using superconducting flux qubits. 



Introduction 

It has recently been demonstrated how an open-ended 
antiferromagnetic Heisenberg spin chain can be used as 
an adiabatic quantum data bus [JJ (see Sec VIA for a 
brief review) . This data bus takes advantage of antifer- 
romagnetic couplings to transfer qubits of information 
adiabatically. First a single qubit, encoded in a single 
spin is joined to an even length Heisenberg spin chain 
slowly enough such that the adiabatic theorem applies. 
Then a single spin on the other end of the chain is sepa- 
rated, again slowly enough for the adiabatic theorem to 
apply. As long as the interactions between the spins on 
the chain are predominately antiferromagnetic, the qubit 
will be successfully transferred from one end of the chain 
to the other. This protocol is illustrated in Fig. [I] An- 
tiferromagnetic spin clusters have been studied for there 
potential usefulness in quantum computing in other con- 
texts, for example in Refs. [21 [3J. 

This paper demonstrates how by applying particular 
unitary operations to spin chains for single qubit gates, 
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Figure 1: Cartoon of an adiabatic quantum bus protocol 
for the Heisenberg spin chain [T|. A spin with the encoded 
qubit is connected to one end of an even length antiferro- 
magnetic chain. Afterwards, the spin on the opposite end is 
removed adiabatically. As long as the chain interactions are 
predominately antiferromagnetic and the adiabatic theorem 
is satisfied the qubit will be transferred. 



and by using a specific spin network for a CNOT gate, one 
can achieve universal holonomic quantum computation. 
This method uses open-loop holonomies, meaning that 
the Hamiltonian is not necessarily returned to the same 
state after the adiabatic process. The methods used here 
can also be extended to non-adiabatic implementations 
of geometrical quantum computing. 

Holonomic quantum computation (HQC) was con- 
ceived and shown to be universal by Zanardi and Rasetti 
j3] and was formulated in terms of a non-abelian Berry 
phase. HQC is considered to be an appealing method for 
achieving fault tolerant quantum computing because of 
its geometrical nature and because it can be implemented 
adiabatically, and therefore has all of the advantages of 
adiabatic quantum computation [5] ■ Although many im- 
plementations of holonomic and geometric quantum com- 
putation are adiabatic, there are examples which are not 

Other proposed architectures for holonomic quantum 
computation use a variety of architectures, including 
superconducting systems with Josephson junctions |S]. 
Further examples propose using quantum dots 01!]. Sin- 
gle molecule magnets have been another system of inter- 
est [7]. A recent proposal has also been made for using 
Holonomies which involve attaching and removing a spin 
from a spin 1 chain [10J. This architecture, although it 
looks superficially very similar to ours, performs com- 
putations locally, at the site of the spin rather than by 
transport as ours does. Ref. [TO] also proposes an im- 
plementation on ultracold polar molecules. Another in- 
teresting proposal involves using a quantum wire with a 
twisted cluster state Hamiltonian [TT]. This proposal is 
similar to ours, but implements the twist in a fundamen- 
tally different way. 

Most other approaches to HQC involve building a sys- 
tem, and explicitly calculating the holonomies caused by 
various manipulations of the system. In our proposal we 
start with a process which has a trivial Berry phasc|26j. 
Real space twists are then performed on the spin chain 
used in this process. Unlike most examples of HQC, all 
results can be derived without explicitly considering the 
curvature of the underlying manifold of states, all single 
qubit gates result from the same underlying Hamiltonian 
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with basis rotations applied to it. 

In principle this idea could be used for more tradi- 
tional, closed-loop holonomic quantum computation by 
twisting a spin chain used to transport the qubit adiabat- 
ically slowly, performing an adiabatic transport protocol, 
untwisting the chain adiabatically slowly, and performing 
the transport protocol in reverse to return to the initial 
Hamiltonian. The net result of this process would be to 
perform a rotation on the initial spin which the qubit 
is encoded and return the spin chain to its initial state. 
In this case spin rotation would be achieved by a non- 
abelian Berry phase acquired via transport on a closed 
adiabatic loop in parameter space, and would therefore 
be a closed-loop holonomic gate. Such a process involves 
unnecessarily many steps however, so it is more conve- 
nient to think of this process in terms of the framework 
of open-loop holonomic quantum computation. 

The mathematical differences of this approach from 
others affords us the advantage that the spectrum of the 
underlying Hamiltonian is the same for all twists, mean- 
ing that, by construction, all single qubit gates can be 
implemented in a way which requires the same anneal- 
ing time to reach a given accuracy. This architecture has 
the advantage that the only operation it ever requires to 
be adiabatically performed is the joining or removal of 
a spin from a chain or cluster. The nature of the twists 
used here also means that a non-adiabatic transport pro- 
tocol could be used instead, and universal computation 
would still be achieved. 

This paper also outlines an implementation of the 
necessary components of this design using supercon- 
ducting flux qubits. Superconducting flux qubits are 
a popular architecture for implementing scalable adia- 
batic quantum computing |12lfl6] . and therefore are a 
natural choice for designing a scalable holonomic quan- 
tum computer. An additional advantage of the use 
of superconducting flux qubits is that the designs tend 
to have spatially extended qubits and a high degree of 
connectivity |1 6| . The large spatial extent of the qubits 
means that a design could be implemented in which a 
qubit would only need to be transferred across a small 
number of spins to be moved from one location in a com- 
puter to any other arbitrary location. For this reason it 
is only necessary that the transport protocol be efficient 
for short chains, as has already been demonstrated in [1], 
rather than in the thermodynamic limit. 

There has been recent experimental work involving 
quantum annealing to degenerate ground state mani- 
folds using currently available superconducting flux qubit 
hardware[17]. In this paper it was demonstrated experi- 
mentally that signatures of quantum behaviors can be ob- 
served in the final state within a degenerate ground state 
manifold. This provides an indication that a ground state 
manifold can be produced accurately enough on current 
hardware that quantum effects dominate over classical 
effects and design inaccuracies. Although the architec- 
ture proposed here cannot be implemented on the hard- 
ware used in |17| , this experiment does provide proof of 



principle for the use of degenerate manifolds in supercon- 
ducting flux qubit systems. 

While it is not the main focus of this paper, we would 
also like to point out that there are other potential meth- 
ods of implementing this architecture. One example of 
such an implementation would be to use a coupled cavi- 
ties scheme similar to the one explored in [18J. For such 
an implementation long spin chains may be required, and 
as a result properties in the thermodynamic limit may be 
important. For such an implementation, the architecture 
given in this paper could easily be generalized to a J1-J2 
spin chain with ^ > 0.25 which is know to be gaped in 
the thermodynamic limit [19 . In this case, another option 
would be to implement the architecture non-adiabatically 
using the methods described in |2Uff2"5] . 

I. SINGLE Q-BIT GATES 

A. The Twisted Spin Chain 

Consider initially an antiferromagnetic Heisenberg spin 
chain. It has been shown that such a chain can act as 
a quantum data bus, both adiabaticallypQ and by us- 
ing the dynamics of its excitations [20-23J. The initial 
Hamiltonian is given by 

JV-l AT-1 

H = = (!) 

i=l i=l 

Now imagine that one inserts a twist into the spin chain 
by applying a local unitary transformation of the form 
x, y, z — > x' , y' , z' on N' = N—L spins[27j, where x' , y' , z' 
are all mutually orthogonal to each other. This yields a 
new Hamiltonian of the form 

L-l AT-1 

Htwist = ^ Oi ■ + <?L • v'l+i + 2-s ®'i ' ^'+1- ( 2 ) 

i=l j=L+l 

Such a twist does not effect the spectrum of the Hamil- 
tonian, and therefore the dynamics of the adiabatic quan- 
tum bus protocol, or other quantum bus protocols which 
may make use of the unitary dynamics of the Hamilto- 
nian. It is important to note, however, that after transfer 
across the chain, the spin will be rotated into the x' , y' , z' 
basis. As we will demonstrate later, transfer through this 
twisted spin chain can perform any desired unitary rota- 
tion on the qubit being transferred, and thus can be used 
to implement any single qubit gate, see Fig. [2] 

One should note that while in this example we consider 
a simple Heisenberg spin chain, gates can be implemented 
in this way on an XYZ spin chain or a J1-J2 spin chain, 
or other sufficiently complex quantum spin Hamiltonians 
|28| . Figuring out which twist to use to perform a given 
gate can be done easily and will be illustrated in the next 
section. 
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Figure 2: (color online) Illustration of a single unitary gate 
implemented by adiabatic transport on a twisted chain. 



B. Example: Implementing a Hadamard Gate 



The local twist to implement the Hadamard gate is 
a x — > <j z , a y — > —errand a z — > a x . This twist can be 
calculated without difficulty, for details see Sec. |VI B| 

One can therefore conclude that the Hadamard gate 
can be implemented by performing the quantum bus pro- 
tocol on the Hamiltonian 



N'-l 

£ 

i=l 



HHadamard = 7 , <7i ' + a AT'°1v"'+l + (3) 



N-l 

a N l<J N' + l - a N' a N' + l + / j a " + 1 



C. Other Single Qubit Gates 

One can perform similar twists to implement any given 
single qubit gate. The calculation to find x' ,y' ,z' in Eq. 
[2]for other gates can be performed in the same way as the 
one in the previous section for the Hadamard. Table ^ 
shows how to implement single qubit gates. These gates 
are sufficient to perform an arbitrary unitary operation 
on a single spin. It is shown in |24| that any unitary ro- 
tation can be approximated to arbitrary precision with 
the gates given in table |T] We have already shown how to 
build an adiabatic quantum bus to move qubit states to 
arbitrary locations in the system. Next we discuss that 
a CNOT gate can be implemented under this architec- 
ture. Then we have demonstrated a universal quantum 
computer. 



II. IMPLEMENTATION OF THE CONTROLLED 
NOT GATE 

A. CNOT design 

Let us now turn our attention to the implementation of 
a controlled NOT (CNOT) using an adiabatic quantum 



Gate Name 


Matrix 


a x ' 




cr 2 ' 


Hadamard 


Ml -1 J 




-a v 


a x 


8 


1 

\ exp(if) J 






a 2 


phase 




-a v 


a x 


a z 


NOT a 




a x 


-a v 


-a 2 



"This gate is needed for the construction of the CNOT 
Table I: Twists for implementing various single qubit gates 



bus protocol. In Fig. |3ji) we show a design for such a 
gate. The time dependent Hamiltonian for this gate is 



Hc7VOt(£; h, tfin) = A(t; tf in )<7 in ■ {a\ + <? 2 + a 3 + a^) 
+d a ia 1 +a 2 +^+a^+h{(al-a z 2 ){l-al) + {al-aX){l+a^)) 
+ (1 - \{t;t fm ))(a' out ■ (at + ff 2 )) + a out (a 3 + a 4 ), (4) 







t < 

o<t<t fvn 



tfin 

1 t> tfi n 



Here "in" refers to the spin which is the input spin, 
where the target qubit is initially encoded; "out" refers 
to the spin to which the target qubit is transferred to, 
"c" refers to the control qubit, "a" to an ancilla to make 
the number of intermediate spins odd. The other 4 spins 
are assigned numbers 1-4. a' out refers to a NOT twist 
being performed on these Pauli matrices, see Tab. [T] 
This gate operates by having 2 channels through which 
a qubit of information can pass. One channel, consisting 
of spins 3 and 4, allows the information to pass through 
the gate unaltered, while another channel, consisting of 
spins 1 and 2 performs a twist on the qubit as it trav- 
els though the gate. The control spin c controls though 
which channel the information travels. The control spin 
is connected with Ising type coupling to spins 1-4 in such 
a way that when the control spin is up the external field 
on spins 1 and 2 cancels with the effect of the Ising bond 
with spin c because (i— < a z c >) = 0, and the infor- 
mation can easily pass though these spins. On the other 
hand (5+ < cr^ >) = 1. So spins 3 and 4 both have 
an effective magnetic field of 2h. For sufficiently large 
h these spins are frozen in the direction of the field and 
will therefore not be able to transport any information. 
As we show in Fig. [3Jd) the net effect is that information 
all travels though spins 1 and 2, and therefore a NOT 
twist is performed. In the case where the spin c is in 
the down direction, information will instead be allowed 
to travel though spins 3 and 4 and blocked on spins 1 
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Symbol 







L7Z 



t 



■1 



Meaning 



spin | control spin 



spin | working spin 



initial AF Heisenberg bond (weak) 



final AF Heisenberg bond (weak) 



fixed AF Heisenberg bond (weak) 



twist to perform NOT operation 



fixed F Ising bond (strong) 



fixed AF Ising bond (strong) 



fixed "up" magnetic field (strong) 



fixed "down" magnetic field (strong) 



Table II: Legend of symbols used in Figs. [3] 



and 2. Therefore in that case the gate acts trivially on 
the qubit. 

Any state of the control spin c can be expressed as 
\ijj c ) = a\ ^) + b\ l) where a and b are complex numbers. 
With an arbitrary input state \ipi n ) = ct\ t) + /3\ I) we 
have before the gate \ip in u) = \ipin) ® \4>c) = ata\ tt 
)+ab| tl)+/9a| lt)+fib\ H). After the gate is performed, 
the final state becomes IV'/m) = o- a \ it) +ab\ f|) +fia\ ft 
) + j3b\ ll) . From these general states we see that 



\4>fin) = 



( 1 o o o\ 

10 
1 

V o o i o J 



\i>t 



(5) 



which is the definition of a controlled NOT gate 



B. Performance of Controlled NOT Gate 

We need to test how this design for a CNOT gate per- 
forms because we cannot rely on previous work to show 
that the qubit is actually transferred accurately. The two 
free parameters in Eq. p4] are the strength of the Ising 
bonds and the fields which we denote by h, and the time 
for the protocol to be performed, tfi n . 

We now examine whether this Hamiltonian actually 
implements a CNOT gate effectively for reasonable val- 
ues of h and t/j n . To test this we need to answer 2 
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Figure 3: (color online) a) Design of a CNOT gate which uses 
the adiabatic data bus protocol. Note that one could replace 
the NOT operation with any other single q-bit unitary, b) 
CNOT system with control central spin up, executes a NOT 
twist on target spin under quantum bus protocol. See Tab. 
[Hi for the meaning of various symbols. Labels are based on 
Eqg] 



questions. Firstly, is the system close enough to the adi- 
abatic limit for reasonable values of tf in 7 Secondly, is 
the desired effect of shutting off one possible path for the 
information achieved for reasonable values of h? To an- 
swer these questions, we examine the overlap of the final 
output state (final state of the "out" qubit in Fig. [3]) 
with the expected output state from a controlled NOT 
gate (Figj4j). Note that because the case where a NOT 
gate is performed, and the case where the gate acts triv- 
ially are related by a simple unitary transformation on 
the Hamiltonian, acceptable performance in one of these 
cases implies acceptable performance in the other. Aver- 
aging over different initial states is therefore unnecessary 
as it would yield the exact same result as any particular 
choice of control and input states. 

Fig. ^demonstrates the effectiveness of this gate. Fig. 
[1] a) shows that for a moderate field and Ising bond 
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III. IMPLEMENTATION USING 
SUPERCONDUCTING FLUX QUBITS 

To build an implementation of this holonomic architec- 
ture based on superconducting flux qubits one needs to 
design circuits that implement Heisenberg spins and the 
appropriate gate couplings between spins. Fortunately 
significant work has already been done, for example in 
Refs. |12f[TB] , towards the design of couplers for super- 
conducting flux qubits. However since the previously dis- 
cussed schemes were based on Ising spin systems, we still 
need to establish a method for designing circuits which 
emulate Heisenberg spins. It is interesting to point out 
that this computational architecture works with the lim- 
ited connectivity of the designs proposed in Refs. |121 - 
ITB"] . Specifically the CNOT gate we have proposed fits in 
a single 2x4 chimera lattice cell like the one used in Ref. 



Figure 4: (color online) Measures of performance of the 
CNOT gate a) one minus fidelity of the output spin versus 
tfi n for h=10 b) one minus fidelity of the output spin ver- 
sus h for tfi„=10 c) output fidelity for initial up spin with a 
NOT performed versus h and tf in light is larger (more posi- 
tive), dark is smaller (more negative) d) gap versus jy— for 
various values of h. The field h and the gap are given in units 
of the Heisenberg bond energy. 



strength the gate can be made to perform well as long 
as the annealing time is sufficient. This figure also shows 
that the gate can continually be made more effective by 
running it longer without having to raise h. The os- 
cillations in the fidelity are related to the time scale of 
small excitations produced during the annealing process. 
In many applications one may have enough control over 
t ji„that the annealing time can be chosen in a way that 
the process lies near one of the local minima of error 
shown in Fig. [5] a). 

Fig. |l]b) shows the effect of h on fidelity for an anneal- 
ing time of 10 (in units of inverse Heisenberg couplings). 
In this figure one can see that increasing h is ineffective 
at improving performance above a certain value. This 
indicates that at this point the field is already effectively 
completely blocking one path that the information trans- 
port can take. Fig. |4]c) shows the combined effects of 
h and i/;„on the output polarization. It is consistent 
with the conclusions we have reached from a) and b) . Fi- 
nally, Fig. |4]d) shows that the system gap remains quite 
large throughout the process. It also demonstrates that 
beyond a certain value of h the gap does not increase sig- 
nificantly with changing h, which is consistent with the 
picture of one path being completely closed to informa- 
tion transfer. It is interesting to note that increasing h 
increases the gap throughout the process. 



A. Flux Qubit Motivation: 

It has been shown in [13] that a qubit which behaves 
like an Ising spin can be constructed from a Hamiltonian 
of the form 



Rising - y^^2C^ +Un 



x\2 



)—U q cos(ai^i) cos(a 2 </>2), 
(6) 



where the applied fluxes <^ act as effective magnetic 
fields, altering the shape of a potential well for the system 
in a way that mimics a spin constrained to move in a 
plane. The constants a simply act to scale the effect of 
the flux. To construct a Heisenberg qubit one needs to 
add a third direction, leading to a Hamiltonian of the 
form 



H 



Heis. 



3 

£< 

71=1 



01 

2C n 



u, 



(7) 



-Ug cos(ai0i) cos(a 2 4> 2 ) cos(a 3 </> 3 ). 

Such a Hamiltonian would allow the shape of the po- 
tential well to be changed along 3 directions and would 
therefore mimic a Heisenberg spin rather than an Ising 
spin. Previously proposed designs also only couple qubits 
along one direction. If a new type of coupler were added 
to the currently implemented circuits which couple the 
spins in the y direction in addition the z direction, then 
an XY model could be implemented. To implement an 
XYZ Heisenberg model, one needs to both design qubits 
which are not constrained to lie in a plane and build 3 
types of couplers, one for each direction in space. 
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IV. FLUX QUBIT DESIGN: 

First consider a CCJJ (Compound-Compound Joseph- 
son Junction) circuit as defined in |13| . 

The effective Hamiltonian of this circuit is given by 



h = £( 



ol 

2C„ 



U n - 



) - U q Peff 



cos 



(8) 

where n G {q, cjj, I, r}. Note that there are more in- 
deces than in jT3| because we do not have the condition 
that cf>L = 4>\ or <t)R = 4> X R . The first two terms of the 
Hamiltonian in Eq(8] are not important for what we are 
trying to demonstrate here [13]. The definition of all of 
these terms can be found in Eq. B4b-f in [TS] , and in Sec. 
|VI C| of this paper. 

Let us make the simplifying assumption that all of the 
critical currents are equal for all junctions. In practice 
there is variability in junction fabrication, but this error 
can be compensated by building a CCCJJ device (see 
Fig. p5j. Let us also assume that our circuit is designed 
in such a way that we can inductively couple the left and 
right loop to each other very strongly such that 4> y = 
b x L = (fif,. These assumptions cause 



f> R and 4>l = 



the equations to simplify greatly (see Sec. VI C ), yielding 



(9) 



13+ = 2/3 L = 2(3 R , 

4ttL 9 / c <f)y 

@HR) = ^ cos (^)- (io) 
This leads to an effective Hamiltonian of the form 



2C n 



K - V n ? 



)-U q /3 + coa(^^-) cos(0 ? ) 

(11) 



When we substitute in /3 + from Eq. 10 this Hamilto- 
nian becomes 



x\2 



tt V^/ Qn j_ T r (<t>n ~ 4»l) 



(12) 



~ Uq 8 ^ 9 Q Ic cos (tp C0 ^Y^ cos (^ z )' 

which is of the form given in Eq. [7] The corresponding 
circuit is shown in Fig. [6] 




Figure 5: (color online) One can build a CCCJJ by replacing 
every Josephson Junction with a pair of parallel junctions in 
the CCJJ. By controlling the flux in any of the smallest loops 
($„ c ) one can effectively change the critical current of the 
junction pair and compensate for manufacturing errors. A 
similar example with a CJJ and a CCJJ can be found in [13] . 




Strong inductive 
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keep f luxes the same 
n smallest 
loops 




Figure 6: (color online) Design using strong inductive cou- 
pling between small loops, where fluxes mimic various mag- 
netic fields applied to a spin. Note that this design assumes 
that all Josephson junctions are identical. 



V. CONCLUSIONS 

We have demonstrated an architecture for a univer- 
sal quantum computer using Heisenberg spin chains and 
clusters. This architecture has the advantage that it can 
be implemented adiabatically and therefore has all of the 
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advantages of adiabatic quantum computing. It has al- 
ready been demonstrated in \Y\ that the single qubit gates 
and data bus used in this computer can be implemented 
with high fidelity for reasonable annealing time. We fur- 
ther demonstrate that a controlled NOT gate can be im- 
plemented at high fidelity with a reasonably short an- 
nealing time and reasonable Hamiltonian parameters. 

In addition to suggesting an architecture, we also pro- 
pose a design to physically realize this architecture. We 
suggest a method for building superconducting flux qubit 
systems which model the low energy degrees of freedom of 
a Heisenberg model. Because the architecture we propose 
can be implemented adiabatically, only the low energy 
degrees of freedom need to be reproduced. We choose 
a superconducting flux qubit implementation because of 
the experimental success of these systems in perform- 
ing non-universal adiabatic quantum computing with an 
Ising spin glass model. Furthermore it has been shown 
that these Ising spin glass models can be realized accu- 
rately enough that the degeneracy of the ground state 
manifolds is not broken. The universal Heisenberg spin 
based computer would represent a significant improve- 
ment over the current non-universal Ising systems be- 
cause it would allow these computers to implement im- 
portant algorithms which the Ising spin glass system has 
not been able to, such as Shors algorithm for factoring 
large numbers [T2HI6] . 
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VI. APPENDIX 




10 ' 1 1 1 1 1 

20 40 60 80 100 

Figure 7: (color online) One minus fidelity of qubit transport 
versus annealing time for quantum bus protocol with linear 
annealing schedule, t is in units of the Heisenberg bond en- 
ergy. 

where A(t, tfi n ) and B(t, tfi n ) are the annealing sched- 
ules, functions chosen such that A(t < 0,t/j„) = 0, 
A{t > t fin ,t fin ) = 1, B(t < 0,t fin ) = 1, and B(t > 
tfi n , tfi n ) — 0- This protocol can also be generalized to 
an XYZ and J1-J2 spin chain (see [T] for more details). 
For simplicity let us consider this protocol with a linear 
annealing schedule such that 

!0 t < 
jf- < t < t fm , (14) 
1 t > tf in 

and 

B(t,t fm ) =1-A(t,t fin ). 

The results of this annealing protocol can be seen in 
Fig. [7] 



A. Review of adiabatic quantum data bus 

The adiabatic data bus protocol discussed here is a pro- 
cess in which an even length anti-ferromagnetic Heisen- 
berg spin chain prepared in its ground state is joined adi- 
abatically slowly to a single spin in a state A spin is 
also removed from the opposite side of the chain, again 
adiabatically slowly. The final removed spin will again 
be in the state The Hamiltonian for this protocol is 



B. Calculation of twist for Hadamard gate 

Showing how to implement any given single spin gate 
using this method is straightforward. Take for example 
the Hadamard gate H, 




(15) 



H A QB{N,t f in) = A(t,tf in )aia 2 + 



N-2 

E 



We now consider the action of T-L on the eigenvectors 
(13) of the Pauli spin matrices, first for a x : 



+B(t,t f 




= z+ 



8 



-> x'_ = Ux. 



Similarly for a z : 



•7o, 



7 = ccjj - (4>° L - 0°,), 



/3_ 7 
7o = - arctan(— tan(-)), 
P+ * 



(17) 



(18) 



(19) 



1 / %/2 I -1 



And for a y : 



P± = Pl±Pr. (20) 
Here we need the additional definitions: 



Ph(R) = Pl(r),+ cos 



2 y pl(r),+ 2 

(21) 



1 [ 1+1 

2 V *— 1 



« + l / 1 



>/2 



^arctan(f^tan(^)), (22) 
Pl(r),+ 2 





In this case the phase factor (</> = f) is irrelevant be- 
cause of the overall U(l) symmetry. 



C. Detailed discussion of simplifying assumptions 
for superconducting flux qubits 



27rLq(Ji( 3 ) ± ^2(4)) 
$0 



(23) 



Let us make the simplifying assumption that all of the 
critical currents are equal, I\ = 1% = I3 = I4. In prac- 
tice there is variability in junction fabrication, but this 
error can be compensated by building a CCCJJ device 
(see Fig. [5|. This assumption causes the equations to 
simplify greatly because Pl{r)- 0, which has the con- 
sequence that <t> L i R \ — > and 7 — > (ftccjj- We can further 
assume that in our design that 4>l = 4>r, this additional 
assumption causes /3_ = and 70 — > 0, which in turn 
causes <jp q — > 0. After this simplification we now have 



The last term U q in Eq. M is a constant which is not 
relevant for this discussion. However the other constants 
in this term are relevant and are defined as follow (Eq. 
B4b-fin[T3]) 



Peff = P+ C0S(-) 1 



'1+C- 



tan(J))2, 



j 8 e // = /8+COB(^-), 



(24) 



(25) 



(16) 
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